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ABSTRACT 

A  study  is  being  made  by  personnel  of  the  Research 
Directorate,  Weapons  Laboratory,  WECOM  to  establish  a 
capability  to  perform  overall  heat  transfer  analysis  for 
any  given  dimensions  of  a  weapon  and  for  any  specified 
propellant  characteristics.  With  respect  to  this  goal,, 
the  propellant  gas  convective  heat  transfer  problem  was 
divided  into  five  subproblems,  (1)  generation  of  thermo¬ 
chemical  properties  for  any  given  propellant;  (2)  transient 
inviscid  compressible  flow  through  the  gun  barrel,  core  flow; 
(3)  transient  viscous  compressible  flow  on  the  bore  surface* 
boundary  layers;  (4)  unsteady  heat  diffusion  through  the  gun 
tube,  and  (5)  unsteady  free  convection  and  radiation  outside 
the  gun  tube  An  exact  numerical  solution  for  the  unsteady* 
compressible,  inviscid  momentum  and  continuity  equations  by 
the  method  of  characteristics  was  obtained.,  In  general*  the 
linear  velocity  gradient  assumption  was  much  more  appropriate 
than  the  uniform  density  assumption.  However,  for  the  first 
portion  of  bullet  travel  (approximately  20  per  cent),  neither 
of  these  assumptions  was  valid  Analytical  boundary  layer 
analysis  procedure  was  developed  on  the  basis  of  the  transient 
compressible  boundary  layer  momentum  integral  equation.  Con* 
vective  heat  transfer  was  evaluated  with  the  Chi  1 ton*Col burn 
analogy  Input,  output,  and  the  computer  program  listing  are 
given  for  the  Lagrange  approximation  and  for  the  exponential 
wall  shear  ' aws  The  free  convection  and  radiation  around 
the  gun  tube  are  of  the  same  order  of  magnitude  The  coupling 
between  the  prob'ems  was  also  discussed 
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I  INTRODUCTION 

As  the  projectile  in  a  weapon  moves  ahead  because  of 
the  high-pressure  gases  created  by  the  burning  propellants, 
the  propellant  gas  will  be  set  into  motion  starting  from 
rest  Since  the  governing  equations  of  fluid  dynamics*  for 
many  problems  of  interest,  are  a  system  of  nonlinear  partial 
differential  equations  dominated  by  real  gas  and  noneauilib- 
rium  effects,  no  general  solutions  exist  in  which  arbitrary 
initial  and  boundary  conditions  are  allowed  Therefore*  the 
flow  field  should  be  examined  and  the  overall  problem  should 
be  subdivided  by  a  consideration  of  the  dominant  features 
only,  i  -e  , ,  (1)  generation  of  thermochemical  properties  for 
hot  gases,  (2)  transient  inviscid  compressible  flow  through 
the  central  section  of  the  barrel,  (3)  unsteady  compressible 
boundary  layer  flow  near  the  bore  surface,  (4)  transient  heat 
diffusion  through  the  multilayered  gun  barrel  walls8  and 
finally  (5)  unsteady  free  convection  and  radiation  outside 
the  gun  tube. 

Constant  gas  (specific  heat  and  gas  constant)  and 
transport  (thermal  conductivity,  viscosity)  properties  are 
common  assumptions  Similarly,  in  the  case  of  problem  (2), 
the  use  of  uniform  density,  uniform  temperature,  and  linear 
velocity  gradient  in  a  gas  flow  between  the  breech  and  the 
bullet  is  popular.  Since  the  governing  equations  of  this 
problem  are  of  the  hyperbolic  type,  they  can  be  solved  by 
the  well-known  method  of  characteri sti cs  Some  of  these 
results  were  presented  at  the  17th  Army  Mathematicians 
Conference,  Redstone  Arsenal,  Ala  None  of  these  assumptions 
are  considered  reasonable  until  after  peak  heating  occurs 
In  general,  however,  the  linear  velocity  gradient  assumption 
is  better  than  the  uniform  density  assumption  Problem  (3) 
can  also  be  interpreted  as  forced  convective  heat  transfer 
In  guns  The  magnitude  of  convective  heat  transfer  in  guns 
is  large  primarily  because  of  the  high  values  of  gas  densities 
that  exist,  and  because  of  the  large  gas=to»wall  temperature 
differences  which  constitute  the  convective  heat  transfer 
driving  potential  Correlation  of  experimental  data  is  quite 
common  with  the  various  analytical  empirical  techniques  The 
most  popular  method  for  heat  transfer  correlation  purposes 
derives  from  the  postulates  of  Nordheim,  Soodak,  and  Nordheim  2 
These  authors  theorize  that  the  propellant  gas  wall  shear 
friction  factor  is  dependent  only  on  gun  surface  roughness 
considerations  The  justification  for  elimination  of  the 
Reynolds  number  as  a  parameter  is  based  on  an  order  of  magni¬ 
tude  estimates  of  boundary  layer  thickness  with  laminar  boundary 
layer  considerations  The  recommended  form  of  the  friction 
factor,  thus,  is  dependent  upon  only  the  ratio  of  surface 
roughness  to  barrel  radius  Consequently,  the  friction  factor 
is  assumed  to  be  Independent  of  space  or  time  within  a  given 
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gun  barrel-  Reynolds'  analogy  is  subsequently  utilized* 
and  the  heat  transfer  coefficient  becomes  simply  propor¬ 
tional  to  gas  density*  velocity,  and  specific  heat-  Total 
calorimeter  data  have  subsequently  been  rationalized  in 
terms  of  the  value  of  the  friction  factor  which  together 
with  internal  ballistic  considerations  obtains  the  spatial 
variation  of  heat  load  to  the  gun  barrel  derived  from  a 
single  shot-  Example  values  of  experimental  friction  factor 
(Cf/2  *  t  /  p  V 2 )  reported  vary  from  0011  to  ,0035t.  Geidt3 

assumes  a  value  of  -002  and  finds  that  his  measured  heat 
flux  is  predicted  generally  within  about  a  factor  of  2, 

Others  have  interpreted  pressure  gradients  within  the  gun 
barrel  in  an  attempt  to  evaluate  wall  shear  for  the  appli¬ 
cation  of  Reynolds'  analogy  for  heat  transfer. 

Another  approach  for  correlation  of  experimental  data 
is  based  on  the  Di ttus -Boel ter  relation  for  steady,  fully 
developed,  turbulent  pipe  flow-  This  relation  is  also  judged 
to  be  inaccurate  in  several  references.  Cornell  Aeronautical 
Laboratory4  proposed  a  combined  analytical -experimental  approach 
based  on  steady  fully  developed  pipe  flow  concepts. 

The  state  of  the  art  in  unsteady  boundary  layers  and 
turbulent  models  is  quite  limited,  Foster-Mi  1 1 er  Associates5 
solved  unsteady  boundary  layers  with  several  approximations  by 
the  integral  method  and  the  method  of  characteristics.  Aero- 
therm6  under  a  contract  with  the  Navy  solved  unsteady  laminar 
boundary  layer  equations  by  the  integral  matrix  procedure, 

Georgia  Tech7  under  a  contract  with  Eglin  Air  Force  Base 
developed  a  solution  procedure  by  the  combination  of  integral 
methods  and  finite  differences  for  turbulent  boundary  layers, 
Aerotherm8  under  a  contract  with  the  Research  Directorate* 

Weapons  Laboratory  at  Rock  Island,  developed  an  analytical 
boundary  layer  procedure  based  on  the  compressible  time- 
dependent  boundary  layer  momentum  integral  equation.  Special 
solutions  of  the  momentum  equation  were  obtained  by  the  classical 
Lagrange  approximation  and  exponential  wall  shear  laws.  The 
solution  was  obtained  by  the  method  of  characteristics.  Con¬ 
vective  heat  transfer  was  evaluated  based  on  the  Chi  1 ton-Col burn 
analogy9  of  the  energy  boundary  layer  to  the  momentum  boundary 
layer. 

Since  the  momentum  and  energy  boundary  layer  equations  are 
not  similar  in  an  accelerating  flow,  an  approximate  solution 
of  the  energy  integral  equation  is  expected  to  yield  better 
heat  transfer  results  than  applications  of  the  Chi Iton-Colburn 
analogy  to  an  approximate  solution  of  the  momentum  integral 
equation.  Another  investigation  is  in  progress  to  determine 
the  validity  of  this  analogy  for  unsteady  boundary  layers  and 
to  obtain  a  modified  form  if  significant  differences  exist  at 
least  for  a  simple  case  such  as  unsteady  incompressible  boundary 
layers. 
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The  unsteady  free  convection  and  radiation  are  discussed 
in  Section  VI.  These  two  contributions  are  of  the  same  order 
of  magnitude.  Since  the  surface  conditions  change  quite 
rapidly  for  automatic  weapons,  the  governing  time-dependent 
partial  differential  equations  are  being  solved  by  numerical 
methods.  The  estimations  based  on  pure  convection  show  that 
the  flow  around  the  gun  tube  is  still  in  the  laminar  range. 

The  solution  of  transient  heat  diffusion  through  gun 
barrel  walls  was  established  quite  satisfactorily  by 
analytical10,  finite  difference11,  and  finite  el ement 1 2  * 1 3  * 1 " 
techniques.  These  are  satisfactory  for  the  purpose  of  estab¬ 
lishing  propellant  gas  convective  heat  transfer  coefficients. 

II.  THERMOCHEMICAL  PROPERTIES  OF  PROPELLANT  GASES 


As  the  propellant  burns,  more  and  more  hot  gases  will 
be  generated.  Typical  gas  in  a  chamber  contains  a  temper¬ 
ature  of  5000°R  and  a  pressure  of  50,000  psi.  The  thermo¬ 
chemistry  of  propellants  involves  determination  of  chemical 
composition  of  gases  (either  by  finite  rate  chemistry  or  by 
chemical  equilibrium  analysis)  and  the  derivation  of  gas 
properties  from  the  composition. 

The  ideal  gas  equation  of  state  does  not  represent  the 
gases  that  exist  in  a  gun  barrel.  For  interior  ballistic 
purposes,  the  id  al  gas  equation  of  state  is  commonly  modified 
by  a  covolume  correction.  The  covolume  term  is  interpreted  as 
the  space  occupied  by  the  gas  when  compressed  to  the  limit. 

This  covolume  term  has  the  empirical  value  of  approximately  1 
cubic  centimeter  per  gram  for  most  of  propellant  gases.  This  is 
insignificant  in  solid  propellant  fueled  rocket  motors,  but,  it 
is  quite  significant  in  gun  interior  ballistic  computations  due 
to  the  high  pressures  involved. 

Covolume  correction  alone  does  not  represent  the  true 
nature  of  the  gas  for  the  entire  ballistic  cycle.  Drastic 
changes  in  chemical  composition,  pressure,  and  temperature 
occurs  throughout  the  firing  cycle. 

Thermochemical  properties  of  typical  propellant  gases 
are  being  predicted  by  use  of  a  NASA-LEWIS  thermochemical 
program  (CEC  71).  The  gases  were  highly  toxic.  (Typical  com¬ 
position  of  the  gases:  Ml 8  (CO  =  0.47,  H 2  =  0.19,  H20  =  0.16, 

N2  =  0.1,  C02  =  0.08)  and  IMR  (CO  =  0.43,  H£  =  0.12,  H20  =  0.21, 

Nj,  =  0.12,  C02  =  0.12)).  IMR  is  better  than  M18  as  far  as  com¬ 
bustion  is  considered,  however,  it  it  still  a  long  way  from  pos¬ 
sible  complete  combustion.  The  consequences  of  incomplete  combust’ 
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are  muzzle  flash,  smoke*  and  fire  in  addition  to  low 
efficiency.  Small  differences  were  observed  between 
expansion  of  frozen  and  chemical  equilibrium  flows.  How¬ 
ever,  quite  significant  differences  exist  for  flame 
temperature,  density,  and  force  of  propellant  between  the 
present  results  (Chemical  Equilibrium  Chemistry)  and  the 
values  reported  in  Reference  15,  which  were  obtained  by 
shortcut  cal culational  methods  Therefore,  future  work 
should  involve  chemical  equilibrium  chemistry  by  consider¬ 
ation  of  all  possible  species  and  reactions  (more  than  one 
hundred);  and,  if  possible,  finite  rate  chemistry  should 
be  included  at  a  later  date 


For  gases  at  moderately  high  pressure,  the  following 
equation  has  been  recommended16? 


PV 

KT  ' 


0,625  b. 


2 

-  + 


0  2869 

- yr 


0,1928  b  •* 

- — yC - - 


where  B(T)  =  Second  virial  coefficient 


and  bQ  =  Vander  Waals  Constant 

Fo*"  powder  gases  at  temperatures  around  2500°C,  the 
Second  virial  Coefficient  is  almost  independent  of  tempera¬ 
ture  A  maximum  value  of  B(T)  (=525  bQ)  is  present  that 

corresponds  to  the  high  temperatures  found  in  powder-gas 
flames  The  value  of  bQ  for  the  powder  gas  mixture  may 

be  estimated  by  the  mixture  rule; 


mi  x 


<»o>i 


where  Xi  -  mole  fractions  of  individual  components 

and  ( bQ ) -  =  Leonard- Jones  bQ  s  for  individual  components, 

1 1 1  CL  AS  S  I_F  I  CAT  I  ON  OF  FLOW  IN  GUN  BARRELS 

The  transport  properties  of  propellant  gases  are  much 
more  important  that  the  thermodynamic  properties  as  far  as 
forced  convection  is  concerned  The  thermodynamic  properties 
and  chemical  composition  of  propellant  gases  were  discussed 
in  Section  II  The  convective  heat  transfer  coefficient  is 
directly  proportional  to  the  thermal  conductivity  of  gases. 
The  thermal  conductivity  of  propellant  gases  that  are  in  use 


cou’d  probably  be  increased  100  per  cent.  Since  the  mean 
free-path  of  gas  molecules  is  inversely  proportional  to 
the  pressure  and  the  number  of  molecules  per  unit  volume 
is  directly  proportional  to  the  pressure*,  the  thermal  con¬ 
ductivity  (of  gases)  that  is  a  product  of  the  two  is  not 
a  strong  function  of  pressure..  However,,  the  thermal  con¬ 
ductivity  is  dependent  strongly  on  the  temperature  of  gases. 
Toward  the  goal  of  establishing  correct  values  the  measure¬ 
ment  of  thermal  conductivity  of  gases  at  elevated  tempera¬ 
tures  is  in  progress  under  a  contract  (DAAF03-71 -C- 0320) 
with  the  University  of  Wisconsin 

Flow  in  a  laminar  boundary  layer  wi'<*  ..veniuslly 
become  unstable  as  the  Reynolds  numbt.  •  •-  .ncreased  The 
boundary  layer  thickness,  skin  friction  and  heat  transfer 
increases  significantly  more  in  turbulent  flow  than  in 
laminar  flow  The  eddy  viscosity  is  the  dominating  mech¬ 
anism  for  such  increase  The  effect  of  eddy  viscosity  may 
be  interpreted  as  if  the  gas  has  a  viscosity  several  orders 
of  magnitude  greater  than  its  molecular  (kinematic)  viscosity. 

The  flow  character i sti cs  are  unknown  for  gun  barrel  flows. 
The  experimental  data  are  ’ack’ng  because  of  the  moving  bullet. 
This  obstacle  may  be  overcome  if  one  takes  advantage  of  the 
similarities  between  the  moving  bullet  (small  mass)  and  a 
moving  shock  in  a  shock  tube  Therefore,  shock  tube  data 
should  be  collected  and  ana’yzed  for  possible  use  in  predicting 
gun  barrel  flow  characteristics 

As  the  propellant  gases  expand  behind  the  projectile,  a 
boundary  layer  forms  at  the  breech  end  and  thickens  as  the 
flow  proceeds  downstream,.  An  unusual  feature  of  the  velocity 
boundary  ’ayer  is  that  it  disappears  as  the  projectile  is 
approached  since  all  fluid  at  the  base  of  the  projectile  must 
be  moving  at  projectile  velocity  Mathemat ica 1 ’y ,  this  amounts 
to  the  requirement  of  an  additional  boundary  condition  at  a 
downstream  location  The  numerical  techniques  applied  to  most 

boundary  layer  problems  cal1  for  'he  specification  of  profiles 
at  the  upstream  end  of  the  flow  i  '  a'low  a  marching"  along 
the  f'ow  direction.  For  the  usual  time-dependent  boundary 
layer  problem*  an  initial  condition  to  describe  the  boundary 
layer  flow  at  time  zero  and  boundary  conditions  as  functions 
of  time  are  required  No  downstream  condition  is  added.  The 
question  logically  arises  then  as  to  how  boundary  conditions 
at  both  ends  of  the  flow  can  be  satisfied  in  any  one  analysis 
Probably,  the  analyses  may  be  carried  out  from  both  ends  of 
the  flow  to  match  the  solutions  somewhere  in  the  middle  of 
flow,  and  thus  the  question  posed  above  may  be  explained  In 
the  past,  three  separate  criteria  have  been  used  to  char¬ 
acterize  fully  developed  turbulent  flow  in  pipes  and  channels. 
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Fully  turbulent  flow  has  been  said  to  exist:  (1)  when  the 
skin-friction  coefficient  is  related  to  the  Reynolds'  number 
by  an  established  law  known  to  hold  at  high  Reynolds'  numbers, 
(2)  when  the  velocity  distribution  in  the  wall  region  follows 
the  well  known  law  of  the  wall,  or  (3)  when  the  flow  is  con¬ 
tinuously  turbulent,  i.e.,  no  intermittence  is  present.  Some 
of  the  results21  indicate  that  these  criteria  yield  different 
Reynolds'  numbers. 

TRANSITION  TO  A  TURBULENT  BOUNDARY  LAYER 


The  laminar  boundary  layer  becomes  unstable  if  the 
Reynolds’  number  becomes  sufficiently  high;  thus  small  dis¬ 
turbances  will  be  amplified  causing  transition  to  a  turbulent 
type  of  boundary  layer. 

For  a  flat  plate  with  zero  pressure  gradient,  the 
laminar  boundary  layer  has  been  shown  experimentally  to  be 
quite  stable  for  length  Reynolds'  numbers  Rex  up  to  about 

80,000  and  can  persist  to  a  Reynolds'  number  of  several 
million  if  the  free  stream  turbulence  is  very  low  and  the 
surface  is  very  smooth.  For  engineering  calculations, 
unless  other  information  is  available,  transition  can 
generally  be  assumed  to  occur  in  the  200,000  to  500,000 
range.  These  figures  are  only  approximate  and  may  be  good 
for  a  smooth  surface  with  a  fair  amount  of  free  stream 
turbulence. 

The  length  Reynolds'  number  Rex  is  inconvenient  for  a 

transition  criterion  since  it  is  based  on  a  constant  free 
stream  velocity  and  may  not  be  meaningful  if  it  is  allowed 
to  vary  with  x  such  as  the  accelerating  flow  in  a  gun  barrel. 

In  such  circumstances,  a  local  parameter  is  preferred  such 
as  momentum  thickness  (<5)  instead  of  x.  If  a  critical 
Reynolds’  number,  Rex,  of  about  360,000  is  chosen  as  a  transi¬ 
tion  criteria  for  a  constant  free  stream  velocity,  the  equi¬ 
valent  criterion  for  accelerating  flow  becomes  Re^  =  0,664 

(600)  »  398.4, 

Note  that  the  Reynolds5  number  based  on  local  distance 
does  not  possess  any  boundary  layer  characteristics,  whereas 
the  Reynolds'  number  based  on  momentun  thickness  does  represent 
the  important  parameter  of  the  boundary  layer.  Since  transition 
to  a  turbulent  boundary  layer  depends  strongly  on  the  growth  of 
the  boundary  layer,  consideration  of  Reynolds'  number  on  the 
basis  of  momentum  thickness  is  logical  to  establish  transition 
criteria. 
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LAMINARI ZATI ON  OF  A  TURBULENT  BOUNDARY  LAYER  IN  GUNS 


The  boundary  layer  is  commonly  assumed  as  turbulent 
for  higher  Reynolds*  number  flows.  Although  the  boundary 
layer  may  be  initially  turbulent,  several  investigations 
of  heat  transfer  from  heated  air  flowing  through  cooled 
axisymmetric  nozzles  revealed  an  apparent  reduction  of 
turbulent  transport  of  heat  along  the  nozzles,  even  though 
the  throat  or  narrowed  part  Reynolds'  numbers  based  on 
diameter  were  as  large  as  4  x  10*  to  2  x  10*.  This 
reduction  In  heat  transfer  is  believ&l  to  be  associated 
with  the  effect  of  flow  acceleration  on  turbulence.  From 
the  low-speed,  essentially  constant  property,  boundary 
layer  measurements  (References  17  and  18)  and  flow 
observations  (Reference  19),  a  turbulent  boundary  layer 
had  become  laminarlike  near  the  wall;  this  was  presumably 
due  to  the  loss  of  turbulent  transport  in  the  wall  vicinity. 
This  process  referred  to  as  laminarization,  apparently  occurs 
when  values  of  a  parameter 


K  * 


due 

”cTx 


exceed  about  2  x  10"*. 

The  effect  of  flow  acceleration  on  the  velocity  and 
temperature  profiles  in  the  nozzle  have  exhibited  the 
following  trends; 

(1)  The  velocity  profiles  become  relatively  flat  in 
the  outer  part  of  the  layer  and  lie  above,  the  1/7  power 
relation, 

(2)  A  laminar  boundary  layer  velocity  profile  fits 
the  measured  profile  satisfactorily  in  the  wall  vicinity. 

(3)  The  temperature  profiles  are  still  in  relatively 
good  agreement  with  the  1/7  power  relation  in  the  outer  part 
of  the  layer. 

(4)  The  thickness  of  the  thermal  layer  is  considerably 
larger  than  that  of  the  velocity  boundary  layer. 

EFFECT  OF  RIFLING 


Separation  of  flow  into  regions  of  low  total -temperature 
near  the  surface  and  high  total -temperature  out  toward  the  free 
stream  by  a  high-velocity  gas  flow  over  a  flat  plate  is  well 
known.  Other  rectilinear  motions  also  demonstrate  such  a 
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separation  of  the  flow  into  regions  of  low  and  high  total 
temperatures.  Such  an  energy  separation  is  much  more 
pronounced  In  vortex  motion  than  in  rectilinear  flows. 

The  super  imposition  of  a  rotational  velocity  on  the 
normal  Poiseuille  flow  creates  lower  total  temperatures  in 
the  core  of  the  flow,  the  effect  becomes  larger  with  increased 
rotation.  Extremely  large  ratios  of  circumferential  to  axial 
velocity  could  create  very  low  total  temperatures  in  the  core 
regions.  The  axial  component  of  velocity  prevents  this 
separation  from  ever  rea<*hing  the  value  obtainable  with  pure 
rotation.  A  greater  separation  effect  exists  in  turbulent 
rotational  flow  than  in  laminar  rotational  flow. 

IV.  UNSTEADY  CORE-FLOW  ANALYSIS 


The  behavior  of  propellant  systems  is  predicted  almost 
entirely  on  the  basis  of  intuition  and  experience,  and  the 
use  of  interior  ballistic  equations.  With  these  equations 
space-mean  averages  are  assumed  for  the  thermodynamic  vari¬ 
ables  and  include  dynamic  effects  only  through  the  inclusion 
of  various  correction  factors  determined  more  or  less, 
empirically, 

In  the  literature,  either  uniform  density  ora  linear 
velocity  gradient  is  assumed  to  be  in  the  flow  field  between 
the  breech  and  the  bullet  at  any  particular  time  during 
weapon  firing  The  objective  of  this  study  is  to  evaluate 
the  validity  of  these  assumptions  in  interior  ballistic 
analyses . 

The  governing  equations,  1.e.9  continuity,  momentum, 
and  energy  of  unsteady,  inviscid,  and  compressible  flow  are 
hyperbolic  in  nature.  In  principle,  one  can  solve  them  by 
the  well-known  method  of  characteri sties .  Along  certain 
(characteristic)  directions,  the  coupled  partial  differential 
equations  can  be  reduced  to  a  set  of  simultaneous  ordinary 
differential  (characteristic)  equations.  The  solution  of  such 
a  set  of  equations  is  further  complicated  by  an  unknown 
nonhomogeneous k  and  unsteady  moving  (bullet)  boundary  con¬ 
dition  and  a  time-dependent  propellant  burning  model. 

The  external  or  wall  friction  and  propellant  burning 
models  are  considered  in  the  formulation  of  the  governing 
equations.  The  Nobel-Abel  equation  of  state  is  used  to 
represent  the  state  of  the  gas.  The  governing  equations  are 
nondimensionalized;  equations  are  derived  for  characteristic 
lines  and  also  for  dependent  variables  along  those  character¬ 
istic  directions. 


8 


As  the  bullet  moves  ahead  because  of  the  high-pressure 
gases,  it  issues  rarefaction  waves  which  will  travel  toward 
the  breech  The  reflected  waves  from  the  breech  will  interact 
with  the  oncoming  rarefaction  waves  and  thus  form  a  network 
similar  to  a  coordinate  grid  system.  The  path  of  these 
waves  may  be  called  characteristic  lines.  The  applicable 
ordinary  differential  equations  along  these  lines  can  be 
solved  by  finite  difference  techniques. 

The  core  flow  model,  while  necessarily  a  gross  over¬ 
simplification  of  the  physical  system,  would  permit  at  least 
a  study  of  the  distribution  of  thermodynamic  variables  and 
provide  more  accurate  interior  ballistics  data  than  the  data 
obtained  by  conventional  methods  to  aid  in  the  design  of 
muzzle  brakes,  noise  and  flash  suppressors,  and  also  gas 
ports  for  automatic  actuating  mechanisms 

MATHEMATICAL  MODEL 

A  model  involving  unsteady,  one-dimensional  motion  would 
apparently  be  the  most  that  could  be  tolerated  In  complexity. 
Moreover,  this  model  does  not  warrant  any  more  complexity  due 
to  the  limited  knowledge  in  solid  propellant  particle  movements 
and  In  associated  burning  models  The  gas  is  assumed  to  be 
inert  and  also  without  viscosity  and  heat  conduction  effects, 
For  convenience,  the  mass  and  energy  are  assumed  to  be  dis¬ 
tributed  continuously  throughout  the  gas  flow,  although  such 
an  assumption  is  unnecessary  to  apply  the  chosen  method  of 
approach  for  the  solution  to  the  problem.  The  gove.iing 
equations  are  given  below 

Conti nui ty 
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Energy ; 


p  y~(e  +  +  pu  ^|-(e  +  +  Q(e  +  *  WQ  *  ^  ^  (3) 


Equation  of  Statg 


p(£  -  n)  *  RT 


(4) 


where 

p  *  gas  density 
u  *  velocity 
A  =  Area 

Q  3  mass  of  gas  produced  per  unit  volume  per  second 
6  *  propellant  density 
p  3  pressure 
f  3  friction  factor 
e  *  internal  energy 
W  =  potential  of  propellant 
n  3  covolume  of  propellant  gas 
R  =  gas  constant 

The  wall  friction  does  not  appear  directly  in  the  energy 
equation  since  its  action  is  that  of  converting  mean  kinetic 
energy  to  thermal  energy 

THERMODYNAMIC  RELATIONS 

Entropy. 


de  3  Cy  dT  *  T  dS  +  pr  dp 


(5} 


or 


dS  .  d£  .  S  +  R 

Z^  p  Z^ 


$■ 


P’(A  •  n> 


(6) 
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(7) 


«p  -  cv  *  Ttlr>p  <It>v 


For  Nobel -Abel  gas,  (yy)  s  (yjp) ^  ^ 


V  i-n 


Cp  -  Cv  *  R  for  Nobel -Abel  gas 


(8) 


Substituting  Equation  8  into  Equation  6  and  Integrating 
the  resulting  equation  yields  the  following  relationship: 


P(J-  -n)Y  *  C  eS/CV 


where 

Cy  =  specific  heat  of  gas  at  constant  volume 
y  =*  specific  heat  ratio 
S  *  entropy 

C  =  integration  constant 
Speed  of  Sound: 

The  square  of  the  speed  of  sound  is  defined  as 


(9) 


(10) 


The  following  equation  can  be  obtained  by  use  of  Equations 
9  and  10: 


rUL 


(11) 


as 


Equation  11  can  also  be  expressed  in  differential  form 


?  _  ]  ~  2pn  d£ 

a  p  1  -  pn  p 


(12) 
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Transformation  of  Governing  Equations: 

The  governing  equations  can  be  conveniently  written 
with  u,  a,  and  S  instead  of  u,  p,  and  e  as  dependent  vari¬ 
ables,  Toward  this  goal,  the  following  equations  are 
derived  by  use  of  Equations  6,  8,  and  12, 


and 


P 


--V-  Pfl  -  [2  4*  .  dS] 

-  1  +  2pn  L  a  Cy 


dp  _  _ 2y  da  1  -  2pn  dS 

p  "  y"-“T  +“2^n  a  ~  y-T ~+  2pn 


(13) 

(14) 


Equation  1  can  be  transformed  into  the  following  form  by  use 
of  Equation  13, 


3a  . 
Fr  + 


u 

3x  T-pn 


a  3u 
7  ’  dx 


a  / 3 S  ,  3S  \ 

7^  {FF  +  u 


7x‘ 


(15) 


*  '  i> 


1  "  1 


Similarly,  the  momentum  equation  can  be  reduced  into  the 
following  form  by  use  of  Equation  14: 


3u  2a  (1  -  pn)  3a  a 
3x  y  -  T  +  2pn  3x  ~  T 


i£&Ill  IS 
TpnT  3  x 


.  «a 


p 


(16) 


The  transformed  energy  equation  is  given  below: 


3S 

3t 


+ 


(17) 


This  is  obtained  by  the  combination  of  Equations  1,  2,  3,  6, 
and  8, 
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Character's  t1_c_Di_recti_oi}s ; 

The  characteri sti c  directions  are  defined  as  the  curves 
along  which  the  derivatives  of  the  fluid  properties  such  as 

ft*  fx*  If*  tx*  If*  lx  are  discont1’nuous°  The  equations  for 

these  curves  can  be  determined  by  consideration  of  the  following 
three  definitions  along  with  Equations  15,  16,  and  17  as  a 
system  of  linear  algebraic  equations  with  variables 


=  du  (18) 


=  da  (19) 


=  dS  (20) 

by  solving  the  above  system 
of  equations  simultaneously  (See  page  14) 

Setting  the  denominator  to  zero  and  performing  algebraic 
manipulations  yields  the  following  equations  for  characterisitc 
di rections 


dx  _  „  .  , 

<Jt  "  u  +  a 

(22) 

3t  =  u  '  a 

(23) 

5?-* 

(24) 

3u  3u  3a  3a  3  S 
cTr'  Ux*  cTt*  ?x*  3t 


x-  and 


3  S 
Tx 


3u,.  .  3Uj v 

37dT  *  Sx*1* 


3a.  .  3a. 

57dT  4  s?** 


as.,  .  SSj. 
TFdT  4  5Sdx 

One  can  obtain  Equation  21  for 


Governing  Equations  Along  Characteristics  Directions: 

The  governing  partial  differential  equations,  namely 
Equations  15,  16,  and  17,  can  be  reduced  to  simple  ordinary 
differential  equations  along  the  characteristic  directions. 
The  reduction  can  be  done  by  setting  the  numerator  of 
Equation  21  to  zero  and  performing  required  algebraic  mani- 
pulatlons  This  procedure  explicitly  uses  the  definition  of 
characteristic  direction  Thus,  the  simplified  governing 
equations  along  the  characteristic  directions  (Equation  22, 
23,  and  24,  respectively)  will  become  equations  25,  26,  and 
27. 
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Solving  for  one  can  obtain  in  determinant  notation: 


(25) 


Du  +  2(1  -  pn)_  Da  =  _i  ,n._P,pJ(.l  -  .2pn),  DS 

157  y  -  1  +  2pn  Dt  C_  y  -  1  +  2pn  cT 


+  §2  (1  -  pn)(W  ♦  4  *  cpT)  +  aQ(P  '  T>  '  “p  - 


Du  2(1  -  pn)  Da 
Ut  -  Y  -T  +  lpn  TTr 


a  11 

7 


* 


-  2pn)  DS 


2prf "  U7 


(1  -  pn)(W  ♦  -  CpT)  -  aQ(-  -  y)  -  ^2  - 


(26) 


Si  .  -2  f„  ♦  Hi  -  e  -  p(l  -  „)] 


(27) 


The  dependent  and  independent  variables  are  nondimen- 
sionalized  with  appropriate  reference  values.  The  new  vari 
ables  are  as  follows: 


u 


A  » 

ao 


_a  r 
ao 


x  7 

l*  T 


aoT  *.  S 

and  S  «  ~  where  L 
L  wn 


length  of  the  barrel 


Propellant  Burning  Model: 

The  propellant  burning  model  is  formulated  by  the 
assumption  that  the  propellant  charge  conforms  to  the  most 
common  burning  rate  equation  in  interior  ballistics,  i.e.. 


r  =  kp 


n 


(28) 


This  is  an  empirical  equation  primarily  based  on  closed 
vessel  tests  such  as  the  bomb  calorimeter.  The  terms  k  and  n 
are  the  ballistic  properties  and  are  available  for  almost  all 
of  the  solid  propellants  that  are  in  use  today* 

The  flow  rate  of  the  burnt  gas  into  the  system  can  be 
obtained  from  Equation  28  by  multiplication  of  it  by  the 
propellant  density,  6  and  the  burning  surface  area,  o. 


r 


I. a,,  Q  =  it? 

(D 

(29) 

T 

Where  to  =  a  /r dt 

0 

(30) 

Boundary  Conditions; 

The  gases  are  at  rest  at  the  breech.  The  gases  at  the 
base  of  the  projectile  are  allowed  to  move  with  the  proj¬ 
ectile  velocity  because  of  a  continuity  requirement. 

As  the  projectile  moves  ahead  because  of  the  high 
pressure  gases  created  by  the  burning  propellant*  the  pro¬ 
pellant  gas  will  be  set  into  motion.  At  every  instant  of 
the  motion*  a  rarefaction  wave  originates  at  the  bottom 
of  the  projectile  and  travels  backward  through  the  compressed 
gas  and  causes  the  gas  to  expand.  However,  from  among  these 
infintely  many  wavelets  only  a  few  are  selected*  to  carry  on 
the  solution  procedure. 

If  the  stepwise  projectile  curve  undergoes  a  change  in 
slope  at  a  point  where  no  incident  wave  occurs,  the  boundary 
conditions  require  that  a  new  wave  be  propagated  into  the  gas 
from  the  point  where  the  projectile  changes  in  velocity. 

The  reflected  wave  from  a  moving  object  may  be  either  a 
rarefaction*  a  compression  or  of  zero  strength*  depending 
upon  whether  the  new  object  speed  is  greater  than*  less  than* 
or  equal  to  the  fluid  velocity  at  the  previous  object  speed. 
Therefore*  the  reflected  wave  will  always  be  a  rarefaction 
wave  in  the  case  of  a  bullet- 

Even  though  continuity  is  required  between  a  bullet  and 
the  gas  at  the  base  of  the  projectile*  this  is  still  an  unknown 
boundary  condition.  However*  a  relation  can  be  formulated  by 
consideration  of  a  dynamic  force  balance  for  the  moving  pro- 
jectile. 


i  »e 


o  £> 


*  PA 


(31) 


where  m  -  mass  of  the  projectile 

v  -  velocity  of  the  projectile 
and  F  =  friction  between  projectile  and  the  bore 

Numerical  Solutions : 

The  interaction  of  rarefaction  waves  originated  from  the 
base  of  the  projectile*  and  the  reflected  waves  from  the  breech 
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are  shown  In  Figure  1.  This  is  constructed  by  the  use  of 
Equations  22  and  23.  in  doing  so,  the  curved  characteri sti cs 
between  the  nodes  are  replaced  by  straight  lines.  This  pro¬ 
cedure  also  yields  the  projectile  path  as  a  function  of  time, 
to  To  accomplish  this  requires  an  iteration  scheme  because 
of  coupling  between  the  governing  equations  and  the  moving 
(unknown)  bullet  boundary  condition.: 

Two  numerical  examples  are  generated,  one  of  which  is 
applicable  to  small  arms  and  the  other  to  an  artillery  weapon. 
These  are  obtained  on  the  assumption  that  the  propellant  is 
burnt  completely  before  the  projectile  starts  moving.  The 
objective  of  introducing  such  a  simplification  is  to  obtain 
the  distribution  of  fluid  properties  along  the  barrel  at  any 
time  with  minimum  effort.  This  assumption  is  unnecessary 
for  the  solution  procedure. 

Only  five  rarefaction  waves  are  initiated  before  the 
reflected  wave  from  the  breech  intersects  with  the  base  of 
the  projectile.  The  number  of  initiated  waves  is  not  a  strong 
function  of  the  net  results;  however,  neither  too  many  nor 
too  few  are  desirable.  A  trial  procedure  is  needed  to  obtain 
more  accurate  results. 

Sufficient  fluid  property  distribution  throughout  the 
barrel  at  any  time  cannot  be  obtained  without  some  inter¬ 
polation  scheme.  Linear  variations  for  velocities  and 
polytropic  relationships  for  pressure,  density,  and  tempera¬ 
tures  are  used  for  interpolation  between  the  two  neighboring 
nodes. 

The  dimensionless  time  versus  the  dimensionless  pressure 
is  plotted  in  Figure  2,  This  solution  is  also  compared  with 
an  approximate  solution  obtained  in  Reference  20.  The  dif¬ 
ferences  are  believed  to  be  primarily  due  to  the  assumption 
of  uniform  pressure  throughout  the  barrel  and  the  use  of  the 
perfect  gas  law.  The  bullet  velocity  versus  distance  along 
the  barrel  is  shown  in  Figure  3,  The  lower  velocities  are 
obtained  as  anticipated  due  to  the  prediction  of  lower  bullet 
base  pressures  than  polytropic  (approximate  solution)  pressures. 
Similar  differences  are  noticed  in  Figure  4  between  the  method 
of  characteristics  solution  and  the  gun  heat  transfer  solution.®* 

The  following  data  on  small  arms  were  used  to  obtain 
Figures  5,  6,  and  7: 

Barrel  diameter  =  1.625  in 

Chamber  length  =  3.0  in 

Charge  mass  =  0,05  lb 
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FIGURE  2  PRESSURE  VERSUS  TIME 


19 


Tl 


TIME (MICROSECONDS) 


2000 


Method  of 

Characteristic  ; 

Gun  Heat  Transfe  • 
Solution 


PROJECTILE  POSITION  L  (FEET) 


FIGURE  4 


PROJECTILE  POSITION  VERSUS  TIME 
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JURE  6  PRESSURE  DISTRIBUTION  IN  SHALL  ARMS 


23 


MICROSECONDS 


VELOCITY  DISTRIBUTION  IN  SMALL  ARMS 


Projectile  mass  =  0.3  lb 

Barrel  length  =  41  In 

Specific  heat  ratio  =  1.2 

Adiabatic  flame  temperature  *  5200°F 
Maximum  pressure  =  33740  psia 

Molecular  weight  =  20 

Obviously,  no  empirical  relation  can  be  devised  for  pressures 
and  density,  at  least  until  the  bullet  leaves  the  muzzle.  The 
data  on  those  curves  represent  time  in  microseconds.  Gas  vel¬ 
ocity  distributions  may  be  approximated  as  linear  after  about 
20  per  cent  of  the  bullet  travel.  Similar  results  on  the 
following  weapon  are  shown  in  Figures  8  and  9: 

Barrel  diameter  s  0.493  ft 

Barrel  length  =  24.97076  ft 

Chamber  length  s  5,56  ft 

Mass  of  the  bullet  *  110,22  lb 

Maximum  pressure  =  13406329.4  psf 

Maximum  density  *  25  lb/ft3 

Covolume  *  0016  ft3/lb 

Specific  heat  ratio  =  1  22222 

In  general,  the  assumption  of  a  linear  velocity  gradient 
is  better  than  the  assumption  of  uniform  density.  Note  that 
the  uniform  density  assumption  implies  a  linear  velocity 
gradient  due  to  the  continuity  equation.  For  the  first  portion 
of  bullet  travel  (approximately  20  per  cent),  neither  of  these 
assumptions  are  good.  For  the  remaining  portion  of  bullet 
travel,  the  linear  velocity  gradient  assumption  is  much  better 
than  the  uniform  density  assumption.  Further  work  is  in  progress 
to  generalize  these  conclusions  by  the  inclusion  of  variable 
propellant  burning  models.  In  this  case,  the  validity  of  these 
assumptions  is  anticipated  to  be  more  questionable  than  in  the 
previous  case. 

V.  UNSTEADY  BOUNDARY  LAYER  ANALYSIS 

The  rate  of  heat  transfer  from  the  hot  propellant  gases 
to  the  barrel  is  controlled  by  the  development  of  the  boundary 
layers.  The  flow  in  the  gun  barrel  boundary  layers  could  be 
laminar,  transitional,  or  turbulent  in  nature.  The  type  of 
boundary  layer  at  a  particular  cross  section  at  any  instant 
need  not  be  the  same  as  at  another  instant.  Since  the  flow 
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has  to  start  from  rest  and  must  also  satisfy  zero  boundary 
layer  thickness  at  the  bullet  base  because  of  the  scraping 
action  of  bullets  laminar  flow  always  exists  in  some  parts 
of  the  gun  barrel  boundary  layers. 

Flow  in  a  laminar  boundary  layer  will  eventually 
become  unstable  as  the  Reynolds'  number  is  increased,  The 
boundary  layer  thickness,  skin  friction,  and  heat  transfer 
'increases  more  rapidly  in  turbulent  flow  than  in  laminar 
flow  The  eddy  viscosity  is  the  dominating  mechanism!  for 
such  increases.  The  boundary  layer  flow  can  be  turbulent 
somewhere  in  the  middle  of  the  flow  between  the  breech  and 
the  bullet  base  A  transitional  regime  should  exist  between 
the  ‘aminar  and  turbulent  regimes.  However,  because  of 
’m’ted  knowledge  about  transitional  regimes,  the  flow  will 
be  assumed  to  change  suddenly  from  laminar  to  turbulent  flow 
at  a  time  and  place  determined  by  a  well-known  lanvinar- 
tu'bulent  transition  criteria.  Therefore,  the  unsteady 
boundary  layer  analysis  is  needed  for  laminar  and  turbulent 
boundary  layers. 

Toward  this  goal,  a  preliminary  study  was  initiated 
with  Aerotherm  Corporation,  Mountain  View,  Californiac 
The  complete  details  of  the  analytical  procedure  developed 
under  this  contract  can  be  found  in  Reference  8.  The 
analytical  procedure  basically  consists  of  an  integral 
techn’que  to  remove  the  dependency  of  the  transverse  in¬ 
dependent  variable  and  the  use  of  the  method  of  charac¬ 
teristics  to  solve  the  resulting  equations.  Simplifications 
such  as  the  Lagrange  approximation,  exponential  wall  shear 
’aws,  constant  wall  temperature,  uniform  density,  pressure 
and  temperature  along  the  barrel,  were  introduced  to  reduce 
the  complexity  of  computations. 

The  energy  equation  of  unsteady  boundary  layers  was 
not  considered  Instead,  the  Chi  1 ton-Col burn  analogy  of 
relating  the  momentum  transfer  to  the  energy  transfer  was 
used  to  obtain  the  convective  heat  transfer  coefficient. 

The  important  equations  of  the  analytical  procedure  given 
above  are  summarized  here  for  convenience. 

Integra1  Momentum  Equation. 


If  *  A.  fj  ♦  (A.  +  A,)e  -  A,  -  0 
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If  the  coefficients  A,  through  A*  can  be  expressed 
as  functions  of  09  x  and  t,  then  a  single  first-order 
equation  for  6  is  obtained.  To  accomplish' thist  in 

the  A«,  term  was  related  to  0  by  use  of  the  following 
turbulence  model : 


..  JS 


=  9  0  <  i  -J*  < 


U;  Jt 


=  5,8 


+  5-2  log  .  y  >  yL 


y^  =  dimensionless  laminar  sublayer  thickness. 

The  results  from  this  model  were  approximated  by  the 
following  relation  which  also  applies  to  laminar  or  turbulent 
condi tions . 


[(1  -  m)^]'  ‘  m 

Reference  properties  were  introduced  to  incorporate 
compressibility.  The  shape  factors  G  and  H  were  evaluated 
based  upon  the  following  assumptions: 


r  -  <£>n 

u  i  o 


T  _  Tw  ,  n  TWxU 

T7  T7  (1 


TT;u  i 


p(~  -  n)  =  RT 


The  quantities  m  and  n  can  be  related  by  the  following 
relation . 


n 


a  2  -  3m 


m 


or  m  3 


2 

nT"J 


The  specific  value  of  G  or  H  then  depends  entirely  on 
T  /T i  ratio.  The  Integral  momentum  equation  along  a 

characteristic  curve 


A, 


can  further  be  reduced  to  the  following  ordinary  differential 
equation: 


+  ( A 2  +  ki)Q  -  Ai,  =  0 

where  D/Ot  represents  substantial  or  material  derivative. 

The  two  equations  may  be  solved  simultaneously  by  ordinary 
methods,  if  the  variations  of  ut  and  p,  along  the  barrel  and 
with  time  are  known.  However,  the  following  assumptions,  in 
addition  to  the  several  mentioned  above,  permit  analytical 
integration  of  those  two  ordinary  differential  equations: 

Longitudinal  gradients  of  gas  and  wall  temperatures  were 
assumed  equal  to  zero 

Longitudinal  gradient  of  gas  density  was  assumed  equal 
to  zero  This  also  implies  a  linear  velocity  gradient 
along  the  barrel  (Lagrange  approximation) 

The  results  of  the  analytical  integration  are  as  follows: 

Location  of  characteristic  lines: 

x  3  xce  "f2(tc)  eMt) 


Momentum  thickness 


-  •  Tx  ft.)’  - 

"  mn[xce  c  3  l - F7TTJ - ^ 


-  m 


where  f, 


:<t>  -  j  Jfdt 
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t  m 


1  -  2m 


r  u  T~rm  f?(t)  rnr 

f 6 ( 1 )  =  J  ”  f4(t)  fg(t)[fi(t)e  2  ]  dt 


J _ 

-  m 


f4(t)  =  f3(t)'  "  EXP  Ij-z-z;  f2(t)3 


f3(t)  =  o  jGVp 


H/G 


m 


f  .  (vl)T-r-5  si 

5  '  v 


L  =  Bullet  Location 
Vp  =  Velocity  of  projectile 

and  prime  represents  evaluation  at  reference  temperature* 
i.e.,  average  of  gas  and  wall  temperature.  Finally,  the 
convective  heat  transfer  coefficient  was  evaluated  by  use 
of  the  Chi lton-Col burn  analogy: 


h  =  pu,  Cp  Pr~2/3  -j 


The  important  quantities  of  interest  are  momentum 
thickness*  Reynolds4  number  based  on  momentum  thickness 
for  laminar- turbulent  transition  criteria,  skin  friction 
coefficient,  and  convective  heat  transfer  coefficient. 

Toe  listing  of  the  computer  program  for  evaluation  of 
these  parameters  among  many  others  is  given  in  Appendix  A. 

The  definition  of  input  and  output  is  available  in 
Appendix  B 

VI.  UNSTEADY  FREE  CONVECTION  AND  RADIATION  OUTSIDE  THE 

guntObT~ 

The  surface  temperatures  of  the  gun  tube  may  reach  as 
high  as  2000°F  at  the  inner  surface  (bore)  and  1000°F  at  the 
outer  surface  because  of  the  high  rate  of  fire.  The  outer 
surface  of  the  gun  tube  is  surrounded  by  ambient  air  which 
may  be  at  a  temperature  between  -60°F  and  160°F  and  may  move 
at  a  veiocity  between  zero  and  60  miles  an  hour.  Since  wind 
velocities  help  to  cool  the  gun  tube  faster  than  free  con¬ 
vection  and  radiation  phenomena,  thermal  design  of  a  gun  tube 
does  not  have  to  Include  wind  velocities.  The  following 
example  illustrates  the  order  of  magnitude  of  free  convection 
and  radiation  effects. 
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Consider  a  unit  length  of  a  gun  tube,  3  inches  in 
outside  diameter,  with  an  outer  surface  temperature  of 
500°F:  Let  the  surrounding  air  be  at  a  pressure  of  one 
atmosphere  and  at  a  temperature  of  70°F.  The  corresponding 
properties  of  air  at  1  atmosphere  and  at  e  film  temperature 
of  285°F  /_  500  +  70*  are  as  follows: 

( - -2—) 

y  *0.08  Ibm/ft-hr 
p  *  0.0533  lbm/ft3 
v  *  0.00041686  f t2/sec 


Cp  *  0.241  BTU/lbm  -  °R 
K  =  0:(.>  BTU/hr-f t  -  °R 

8  =  T~  =  7XS  rf 


'film 


Grashof  Number  =  Gr  =  DiP.ffi.AT  *  1.66956  X  10s 

yC 

Prandtl  number  =  Pr  ■  *  0.964 

Gr.Pr  *  1.60946  X  10* 

This  dimensionless  parameter  lies  in  the  range  of  10“  to  10*. 
By  use  of  the  approximate  formula  given  by  McAdams,  the 
Nusselt  number  can  be  calculated  as 

Nu  =  0.59  (GroPr)1/"  =  21-015 

Then,  the  heat  transfer  coefficient  becomes 

h  *  Nujj  *  1.68117  BTU/hr-f t-°F 

The  rate  of  heat  loss  per  unit  length  of  the  gun  tube  becomes 


Q  * 
Mc 


h  (ttD) AT 


*  567.7668  BTU/hr-ft 

The  rate  of  energy  emission  from  the  gun  tube  to  the 
surrounding  air  is  given  by 

Qr  *  eoAFJ2Tw- 
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where  e  -  Emissivity  of  the  surface 

o  *  Stefan-Bolt2mann  Constant 

A  s  Emission  surface  area 

F  j  -  Radiation  interchange  factor 

Tw  *  Gun  tube  outer  surface  temperature 

The  emissivity  of  the  surface  is  assumed  to  be  unity. 
This  is  desirable  to  maintain  low  temperatures  in  a  gun 
tube.  Perhaps*  this  could  be  achieved  by  an  oxidized 
coating.  Since  the  surface  is  convex*  the  interchange 
factor  F . 2  is  assumed  to  be  unity.  Also*  energy  absorption 
is  assumed  to  be  zero  because  of  the  absence  of  other 
radiating  bodies. 

Therefore,  the  net  heat  loss  by  radiation  per  unit 
length  of  the  gun  tube  becomes 


Qr  *  oA(V  -  V) 

*  0.171  X  1 0'  * (itD)  ( 960*’  -  530" ) 

-  1039  BTU/hr-f t 

The  boundary  layers  which  arise  by  free  convection 
a>so  change  to  turbulent  flow  where  they  have  reached  a 
certain  thickness  In  air.  the  change  occurs  at  a  critical 
Grashof  number  around  Grx  1  10s,  This  corresponds  to  a 

ue6 

Reynolds"  number  Refi  =  ■  550.  The  numerical  example 

cited  above  is  still  ’n  the  laminar  range, 

VII  COMBIhEO^ANALYSIS 

The  unsteady  heat  transfer  analysis  for  the  chosen 
ammunition  and  gun  was  divided  into  five  problems  in 
Section  I,  The  solution  to  these  problems  was  discussed 
in  Sections  II  through  VI.  The  combined  analysis  to  ful¬ 
fill  the  overall  objective  is  now  in  order.  The  solutions 
obtained  for  the  individual  fictitious  problems  do  not 
represent  the  solutions  for  the  real  gun  tube  problems 
because  of  continuous  change  in  boundary  conditions.  A 
method  of  predicting  and  correcting  these  boundary  con¬ 
ditions  is  desirable  Compatible  boundary  conditions  must 
be  introduced  at  interfaces  between  the  problems- 
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Since  Problem  10  mentioned  in  Section  I,  is  a  huge 
one,  it  was  decided  to  generate  ahead  the  thermochemical 
properties  for  various  pressures  for  any  chosen  propellant. 
This  information  will  be  fed  in  as  input  in  the  form  of  a 
table  for  Problem  2.  An  interpolation  scheme  will  be 
introduced  to  select  the  appropriate  quantities. 

To  solve  Problem  2,  one  has  to  know  the  location  of 
boundaries,  in  particular,  the  Interface  between  Probelms 
2  and  3„  Since  the  location  of  this  interface  is  unknown 
and  also  the  thickness  of  boundary  layers,  in  general,  is 
small,  one  can  reasonably  assume  that  the  interface  is  the 
bore  surface  for  Problem  2.  However,  one  should  include 
the  information  at  the  interface  such  as  the  heat  loss  to 
the  gun  tube  and  skin  friction  into  the  governing  equations. 
These  quantities  will  be  dependent  upon  the  development  of 
boundary  layers,  i.e..  Problem  3  and  also  the  temperature 
of  the  bore  surface.  The  outer  edge  of  the  boundary  is 
unknown.  Sometimes,  the  assumption  that  the  outer  edge 
extends  to  the  centerline  of  the  tube  (without  any  loss 
of  generality)  is  convenient. 

Even  though  the  boundaries  are  well  defined  for  a 
chosen  gun  for  Problem  4,  the  information  at  the  boundaries 
such  as  heat-in  due  to  forced  convection  and  heat-out  due 
to  free  convection  and  radiation  (Problem  5)  are  lacking 
because  these  in  turn  are  dependent  upon  Problem  4. 

It  Is  not  clear  at  this  time  how  strong  the  inter¬ 
action  or  coupling  is  between  the  problems,  However,  it 
is  anticipated  that  for  any  reasonable  time  increment,  one 
can  start  with  Probit.?  2  and  proceed  to  the  end  of  Problem 
5  by  solving  one  after  She  other  and  using  as  much  of  the 
latest  information  as  possible.  Of  course,  one  has  to  use 
the  interface  boundary  information  from  the  preceding 
problem  This  procedure  may  be  satisfactory  if  no  strong 
interaction  exists  between  the  problems  and  if  small  time- 
increments  are  chosen.  Otherwise,  one  has  to  iterate  the 
procedure  given  above  until  the  results  converge.  The 
possibility  of  decoupling  Problem  2  and  5  from  the  procedure 
above  should  be  investigated 

Cookoff  is  still  a  problem  with  caseless  ammunition 
and,  to  some  extent,  with  conventional  ammunition,  if  the 
chamber  is  designed  for  reduced  weight.  Accurate  deter¬ 
mination  of  chamber  temperatures  is  essential  to  satisfy 
the  desired  reduced  weight  and,  at  the  same  time,  to 
eliminate  cookoff  hazards. 
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VIII 


CONCLUSIONS 


A  study  was  initiated  to  establish  a  capability  to 
perform  overall  heat  transfer  analysis  for  any  given 
dimensions  of  a  weapon  and  for  any  specified  propellant 
characteristics,  Toward  this  goal,  the  propellant  gas 
convective  heat  transfer  problem  was  divided  into  five 
subproblems;  (1)  generation  of  thermochemical  prop¬ 
erties  for  any  given  propellant*  (2)  transient  inviscid 
compressible  flow  through  the  gun  barrel  (core  flow), 

(3)  transient  viscous  compressible  flow  on  the  bore 
surface  (boundary  layers),  (4)  unsteady  heat  diffusion 
through  the  gun  tube,  (5)  unsteady  free  convection  and 
radiation  outside  the  gun  tube 

The  prediction  of  the  composition  of  the  propellant 
gases  by  chemical  equilibrium  chemistry  for  typical  small 
arms  propellants  such  as  M18  and  IMR  revealed  that  the 
gases  (more  than  40  per  cent  carbon  monoxide)  were  highly 
toxic.  The  flame  temperature,  density,  and  force  of 
propellant  differ  significantly  from  the  results  reported 
in  AMCP- 706 - 1 50 

The  classification  of  flow  in  gun  barrels  is  as 
important  as  the  analysis  on  any  one  of  the  problems 
mentioned  above  The  boundary  layer  flow  is  turbulent 
for  most  of  the  region  inside  a  gun  barrel  The  chances 
for  relaminar wation  are  remote  The  rifling  may  Increase 
the  heat  transfer  as  much  as  25  per  cent 

An  exact  numerical  solution  for  the  unsteady,  com¬ 
pressible.  inviscid  momentum  and  continuity  equations  by 
the  method  of  characteristics  was  obtained  In  general, 
a  linear  velocity  gradient  assumption  is  found  to  be  much 
better  than  the  uniform  density  assumption  For  the  first 
portion  of  bullet  travel  (approximately  20  per  cent), 
neither  of  these  assumptions  are  satisfactory  Further 
work  is  in  progress  to  incorporate  the  burning  model  and 
energy  equation  into  the  governing  equations,  and  to 
obtain  the  complete  solution  by  the  method  of  characteristics 

Analytical  boundary  layer  analysis  procedure  was 
developed  based  on  the  transient  compressible  boundary 
layer  momentum  integral  equation.  Convective  heat  transfer 
was  evaluated  by  use  of  the  Chi  1 ton-Col burn  analogy.  The 
Input,  output,  and  listing  of  computer  program  ware  given 
for  the  case  of  Lagrange  approximation  and  exponential  wall 
shear  laws.  An  investigation  to  verify  the  validity  of  the 
Chi  1 ton-Colburn  analogy  (relating  momentum  transfer  to  energy 
transfer)  of  steady  flow  to  unsteady  flow  is  in  progress. 
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The  analysis  of  transient  heat  diffusion  through 
single  or  multilayered  gun  tubes  was  developed  under  an 
In-House  Laboratory  Independent  Research  project,,  This 
Is  satisfactory  for  the  present  i nvesti gation.  The 
amount  of  free  convection  and  radiation  is  found  to  be 
of  the  same  order  of  magnitude.  The  analysis  of  forced 
convection  on  the  outer  surface  of  the  gun  tubes  because 
of  severe  wind  conditions,  is  unnecessary  due  to  the 
favorable  resulting  effects.  The  boundary  layers  of 
typical  free  convective  flow  on  the  outer  surface  of  the 
gun  tube  are  estimated  as  being  of  the  laminar  type.  The 
analysis  of  the  free  convection  and  the  radiation  problems 
by  explicit  finite  difference  methods  are  in  progress.  The 
interaction  between  the  problems  and  the  overall  procedure 
to  combine  all  these  problems  was  discussed. 
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APPENDIX  A 

Listing  of  Computer  Progrem 


//MGHU1CH  JOB  ( YYY YYY YY , 5 , bCCC > , *  PUNCH  THAT  DATA* 

C 

C  THIS  IS  THE  MAIN,  JR  CALLING  PROGRAM  FOR  THE  SOLUTION  UF  THE  TIME 
C  DEPENDENT,  COMPRESSIBLE,  INTEGRAL  BOUNDARY  LAYER  MOMENTUM  EOUAT ILN 
C  BY  THE  METHOD  UF  CHARACTER  I STICS.  THE  DEVELOPMENT  OF  THE  METHOD,  AND 
C  THE  EQUATIONS  EMPLOYED  FUR  THIS  SPECIAL  SOLUTION  APPLICABLE  TO  GUN 
C  BOUNDARY  LAYERS  ARE  PRESENTED  IN  AEROTHERM  REPORT  70-18  -  T.J.  DAHM 
C 

COMMON  I ,  I  I 

COMMON  DOTS2GO<,F1L200<*F2%200<*F4S200<»F6S20C<»REFRS200<»RNLZ200< 
2,SFHS2GC<»TS20G<,XFT*2C0<, VK*20C<» T1 *200<,XI C220< 

COMMON  AC, ALFA, BETA, DTUB, F  INTS,FIS,PRDL,RO,  SFG,TW,  V, VKX  ,  VMX 1 ,  VM 

2X2,VMX3,VMX4,VMX5,VNU,DT,ELD 
OMEGA#. 76 
C 

C  READ  BASIC  INPUT  DATA 
C  DTUB  #  GUN  TUBE  DIAMETER  IN  INCHES 

C  ELZ  #  CHAMBER  LENGTH  MEASURED  FROM  THE  BREECH,  IN  INCHES 

C  ELD  #  LOCATION  OF  DISCONTINUITY  FROM  THE  BREECH,  IN  INCHES 

C  TW  #  WALL  TEMPERATURE,  IN  DEG.  R 

C  CP  «  SPECIFIC  HEAT  AT  CONSTANT  PRESSURE, IN  BTU/LBM/OEG. R 

C  II  #  NUMBER  OF  PROJECTILE  PUSITIONS  TO  BE  CONSIDERED 

C  PRDL=PRANDTL  NUMBER 
C 

READS  5 , 1<  DTUB, ELZ, ELD, TW, PR DL, CP, II 

1  FORMA TS6E 10.3, II 0< 

ELDtfELD/ 12. 

PRFC#t • 7 1/PRDL<** .666  7 

READ  AND  PROCESS  MURE  INPUT,  AND  INITIALIZE  INTEGRATIONS 

CALL  IN  1 T 
DO  21  I«l,Ii 

READ  INTERNAL  BALLISTIC  DATA  IN  CHRONOLOGICAL  ORDER 
T  *  TIME  IN  SEC. 

RO  #  GAS  DENSITY  IN  LHM/CU.F T . 

Tl  H  GAS  TEMPERATURE,  IN  DEG.  R 

X  #  PROJECTILE  POSITION  RELATIVE  TO  INITIAL  POSITION,  IN  INCHES 

V  #  PROJECTILE  VELOCITY,  IN  FT./SFC 

RE  ADZ  5,2<T«I<,R0,T  UI<,X,  V 

2  F0RMATZ5E IG. 3< 

TZI<#TSI</1000. 

XFTS1<#ZXCELZ</I2. 

IFSI-K  10,10,11 

10  0T#0, 

GO  TO  12 

11  0T#T*I<-T«I-1< 

CALCULATE  VISCOSITY  ZLBF/HR/F  T<»  THERMAL  CONDUCTIVITY 
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C  *BTU/HR/FT/DEG.R<,  KINEMATIC  VISCOSITY,  AND  REYNOLDS  NUMBER 

12  VMU#.044**T1*I</S30.<**0MEGA 
VK«I<#CP*VMU/PRDL 
VNU4VMU/3600./R0 
RN#V*DTU8*RO/VMU*3CiO. 

RNL*I<#RN*X/DTUH 

FORM  CHARACTERISTIC  FUNCTIONS 

CALL  CHARF 

IFIX.EQ.C.OIGO  TO  211 

HLAM  #  PROJECT 1LE  BASE  HEAT  TRANSFER  COEFFICIENT  -FROM  EQUATION  113 
OF  AEROIHERM  REPORT  70-lfi,  IN  BTU/SEC/  SQFT/DEG.R 

HLAM#.0022*PRFC*CP*VMU/X*RNL  *1 <* 12.  /3600./.796 
GO  TO  212 

211  HLAK=0.0 

2 12  CONTINUE 

PRINTOUT  TRAJECTORY  AND  OTHER  DATA-  SAME  UNITS  AS  INPUT 
WRITE*  6,3<  T*1<»RU»T1*I<*X,V»HLAM 
3  FORMAT* 10X,  6*1X,E12.S« 

21  CONTINUE 

OBTAIN  CHARACTERISTIC  SOLUTIONS 

CALL  GUNBL 
STOP 
END 

SUBROUTINE  INIT 
COMMON  I, II 

COMMGN  0UTt20G<tFl*200<»F2*200<»F4*200<,F6*20G<,REFR*20G<»t<NL*200< 
2»SFH*200<,T*2C0<,XFT*2C0<,VK*2CC<,T1*200<,XIC*20< 

COMMON  AC,ALFA,BETA,DTU8,F INTS ,F 1 S , PRDL , RO , SFG »TW,  V,VMX ,VMX l,VM 

2X2»VMX3»VMX4»VMXS»VNU,DT,fcL0 
9  F0RMAT*3E1G.3< 

READ  INPUT  PAR4MFTERS  M*#VMX<  AND  A*#AC<  OF  EQUATION  SI,  REPORT  70-18 
AND  SHAPE  FACTOR  G.  THESE  ARE  ALL  DIMENSIONLESS 

READ*S,9<VMX,SFG,AC 

FORM  EXPONENTS  THAT  WILL  BE  USED  LATER 

VMXl#l.-VMX 
VMX2*2.*VMX 
VMX  VMX/VMX 1 
VMX4IU1.-2.*VMX</VMX1 
VMXS# 1 . / VMX 1 
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INI  IALIZE  INTEGRATIONS 

F1S#0. 

FINTS#0. 

F1Z1<#0. 

F2Z1<#0. 

F4<1<#0. 

F6ZK/U0. 

RETURN 

END 


SUBROUTINE  CHARF 

THIS  SUBROUTINE  FORMS  THE  CHARACTERISTIC  FUNCTIONS  F1,F2,F3,F4,  AND 
F6i  EQUATIONS  B2, 81, 8C, 79,  AND  122  RESPECTIVELY  OF  REPORT  70-18 

COMMON  I,  II 

COMMON  DOT*200<,F1*200<.F2*200<,F4*200<,F6*200<,REFR«200<,RNLT;200< 
2fSFH^200<»TZ200<tXFT*200<fVKZ2CO<»TlZ200<»XICZ20< 

COMMON  AC, ALFA, BETA, DTUB, PINTS ,F IS, PRDL,R0, SFG,TW,  V , VMX , VMX I , VM 

2X2 , VMX3, VMX4 »VMXb,VNU»DT ,FLD 
11  FORMATS 11Z1XE11.4<< 

OBTAIN  BOUNDARY  LAYER  SHAPE  FACTORS 

CALL  COMBL 
IFZI-K  34,  34,  33 

33  F1«I<#V/XFTZI< 

F2*I<#F2*I-1<&SFUI</SFG£F1S</2.*DT 

EXF#SFHZI</SFG 

F3#RU*SFG*V**EXF 

F4iI<#F3**VMX5*EXP<2./VMXl*F2£l« 

PINT  #F4Z I<*VNU**VMX3/ SFG*?F 1ZI<*EXPZF2ZI «<**VMX4*RFFRSI  < 
F6XI<#F6ZI-1<CZFINTGFINTS</2»*DT 
F1S#FIXI</SFG 
F INT  S#F INT 

PRINTOUT  INTERMEDIATE  DATA  FOR  RFFERENCF  PURPOSES 

34  WRITE'  6,  ll<T«I<»FrU<,F2<I<,F4*I<,F6*I<,RNHI<,XFTtI<,D<Tm<,SP 
2H*I<»SFG»REFR-SI< 

RETURN 

END 
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SUBROUTINE  COMBL 


THIS  SUBROUTINE  EVALUATES  THE  BOUNDARY  LAYER  SHAPE  FACTORS  OF 
FIGURE  9  OF  REPORT  70-18,  VALID  ONLY  FOR  1/7  POWER  AND  TEMPERATURE 
PROFILES.  EVALUATIONS  BASED  ON  CURVE  FITS  TO  EXACT  RESULTS. 

COMMON  I, II 

COMMON  DUTX200<*  F1?200<,F  2£200< *  F4X200<»  F64200<»R£FR:?200< , RNLX200< 
2,  SFH3>200<»  T?2C0<»  XFT720G<»VKt200<»  T1 ?200<»XI CS20< 

COMMON  AC, ALFA, BETA, OTUB, F 1NTS,F IS, PRDL , RO» SFG,TW ,  V, VMX , VMX 1 , VM 

2X2,VMX3,VMX4,VMX5,VNU,DT,ELD 
CALL  PROPS 

SFGrfl.2857C.C545*BETA-.0369*BETA**1.17 
SIGMA# 1 . -ALFA 

SFHSI<#1.28S7-1.25  72*SIGMA-.025*SIGMA**2.2I> 
D0TU<#ALFA**10.2857C8.06857*BETA< 

RETURN 

END 

SUBROUTINE  PROPS 

THIS  SUBROUTINE  EVALUATES  THE  BOUNDARY  LAYER  REFERENCE  PROPERTY 
DATA. 

COMMON  1,11 

CUMMCi.N  t)OT*20(J<,FU20L<,F?%200<,F4*200<,F6*200<,RtFRT:20C<,RNL,C2CO< 
2,SFH«2G0<,T-6200<,XFT*?CG<,VKS200<,T1*200<,XIC£2C< 

COMMON  AC, ALFA, BETA, DT UB , F I NTS,F 1 S , PRDL , RO, SFG, T W ,  V, VMX , VMX 1 , VM 

2X2 ,  VMX3,  VMX4 ,  VMXt>,  VNU,DT,ELD 


COV#GAS  CUVOLUMESCUFT./LBM<.  IT  IS  EQUAL  TO  0.  FOR  THIS  SAMPLE  CASE 

COV#C. 

OMEG#.7A 

ROCO#RO*COV 

TR#TW/T1*1< 

TRP#1.-TR 

ALFA  AND  BETA  ARE  PARAMETERS  NEEDED  TO  E VALUATF  BOUNDARY  LAYER 
SHAPE  FACTORS 

ALFA#TRCROCn*TRP 

beta#*trp**i.-roco«/alfa 

TPmW£TUI«/2. 

TR#TP/TUI< 

ROR# I ./ vROCOEHl .-ROCO<*TR< 

VMUR#TR**OMEG 
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REFR  IS  EQUATION  B40F  REPORT  70-18 


REFRS I<#R0R**VMX4*VMUR**VMX3 
RETURN 
END 

SUBROUTINE  GUNBL 

THIS  SUBROUTINE  OBTAINS  THE  CHARACTERISTIC  SOLUTIONS  AFTER  ALL  IN¬ 
TERNAL  BALLISTIC  DATA  HAVE  BEEN  PROCESSED 

COMMON  I, II 

COMMON  DOT*200<,FU200<,F2i200<»F4*200<,F6*200<»REFR«20C<tRNLX200< 
2»SFHS200<»  TS200<»XFTS200<»VKS2C0<#Tl'S200<»XIC3!2O< 

COMMON  AC,ALFA,BETA,DTUR,FINTS,F1S,PRDL,R0,SFG,TW,  V, VMX , VMX 1 , VM 
2X2,VMX3,VMX4,VMX5,VNU,0T,ELb 
1  FORMATS I  1,  E 12. 5< 

6  FORMATS  1 1C< 

7  F0RMATS8E10.  3< 

8  F0RMATS4S1XF7.3C, 1XF 7.0, 1XF7.2, 3SF8.3C, 1XF7 . 2 , IXF7 . 5 ,2S1XF9. S<, 2 41 
2XF7.4<< 

9  FORMAT  S1H1.44X5HTIME#,  F7«,3»  1X12HMILLI SEC0NJS/30X 
220HPR0JECT ILE  PUSITI0N#,F?.3,lX6HINCHES//< 

10  FORMATS  4X2HTC6X2HXC7X1HX6X3HHTC4X4H0D0T3X8H0/SQ-DB<?X5HTHETA 
23X5HDELT  A2X6HDELSTR3X3HSHR4X4HSHRT4X4HTHTT4X5HRTHTT4X4HNUTL5X2HX I < 

READ  INPUT  FOR  SOLUTION  OUTPUT  TIMES  AND  STATIONS 

NXT#0  IF  DESIRED  PRINTOUT  POSITIONS  HAVE  NOT  BEEN  DEFINED 
NXT#  1  IF  DESIRED  PRINTOUT  POSITIONS  HAVE  BEEN  DEFINED 
NXT  GREATER  THAN  1  WILL  CAUSE  RETURN  TO  THE  MAIN  PROGRAM 
TIME  IS  THE  SOLUTION  TIME,  IN  SECONDS 

15  READI5, IINXT.TIME 
IFSNXT-1<20,22,26 
20  FPRT#SAC/VMXK**VMX5 
C 

C  THE  FOLLOWING  INPUT  OFFINES  PRINTOUT  POSITIONS.  NN  IS  THE  NUMBER 
C  OF  POSITONS  BETWEEN  THE  PROJECTILE  AND  BREECH,  AND  XIC  SBETWEEN  0. 

C  AND  l.<  DEFINES  THE  DIMENSIONLESS  ORIGIN  IN  SPACE  OF  CHARACTERISTICS 
C  ALONG  WHICH  THE  OUTPUT  DATA  ARE  EVALUATED.  THESE  DATA  ARF  READ  ONLY 
C  IF  NXT #0 
C 

RE A0( 5, 6  INN 

READS  5. 7<SXICSN<»N#1,NN< 
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22  DO  24  i#l»  1 1 


C 

C  FIND  CHARACTERISTIC  AND  OTHER  RFQLIRED  FUNCTIONS  AT  THE  OUTPUT  TIME, 

C  BY  INTERPOLATION 
C 

DELT#TIME-TSI< 

IF XDELT<23, 23,24 

23  FRCC#-DELT/£T1!I<-T'«I-1<< 

FRCH.-FRCC 

F1T#FUI-1<*FRCCCF1*I<*FRC 

F2T#F2*I-l<*FRCC£F2iI<*FRC 

F4T*F4XI-1<*FRCCSF4*I<*FRC 

F6T#F6*I-1<*FRCCSF6?I<*FRC 

RNLT  #RNL4I~1<*FRCC£RNL?I<*FRC 

XFTT#XFT4I-1<*FRCC£XFT?I<*FRC 

DOTT  #DOTtI-l<*FRCC£OOT?I<*FRC 

SFHT#SFH2I-1<*FRCC&SFH*I<*FRC 

REFRT#REFR< I- i<*FRCC£REFRZ I <*FRC 

TlT#TltI-l<*FRCCCTl'tI<*FRC 

VKT#VKXI-1<*FRCC£VK*I<*FRC 

XFTIN#XFTT*12. 

TMS<KTIMC*1000. 

C 

C  PRINT  OUTPUT  TIME  AND  PROJECTILE  POSITION  AT  THAT  TIME 
C 

WKITEI6,  9)TMS,XFTIN 
WRITEI6, 10) 

RNLM#RNLT**VMX 

XOD#XFTIN/DTUB 

DBF#.023**RNLT/X00<**.B*PRDL**.4 
GNFtfRNLT**VMXl*PROL**.3333/xaO 
GO  TO  25 

24  CONTINUE 

25  OU  35  N#1,NN 

f 

V* 

C  FIND  ^BETWEEN  HERE  AMO  STATEMENT  3i<  THE  CHARACTERISTIC  ORIGIN  IN 
C  SPACE  *XC<,  AND  THEN  THF  CHARACTERISTIC  ORIGIN  IN  TIME  *TC<,  BASED 
C  CN  THE  DIMENSIONLESS  POSH  ION  ORIGIN  OF  THE  OUTPUT  CHARACTERISTIC,  XlC 
C 

XC#XFTT*<1.-XIC.SN« 

C 

C  CHECK  FOR  CHARACTERISTIC  ORIGIN  UPSTREAM  OF  THE  DISCONTINUITY 
C 

DX#XC-£LD 
i F?DX<  30,30,21 
C 

C  DETERMINE  iF  CHARACTERISTIC  ORIGINATES  IN  THF  CHAMBER 
C 

21  0X0XC-XFTS1< 

I FXUX<  29,29,31 
29  TC#0. 

F2TC#0. 

F6TCOO. 

XFAC#XC 
GO  TO  34 
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STATEMENT  30  THROUGH  28  DOES  SPECIAL  CALCULATION  TO  CONSIDER  A  DIS¬ 
CONTINUITY  IN  THE  GUN  BARREL  ACCORDING  TO  BOUNDARY  CONDITION  125  OF 
REPORT  70-18 

30  TC#*ELD-XC</EL0*TIME 


XCrfELD 

DO  28  I#1,II 

DELT#TC-T*I< 

IF*DELT<  27,27,28 

27  FRCC#-UELT/*T*I<-T*I-1« 

FRC# i.-FRCC 

F2TCMF2* I- l<*rRCC£F23I <*FRC 
F6TC#F6?I-1<*FRCC£1F6*I<*FPC 
XFAC#XC/EXP*F2TC< 

GO  TO  34 

28  CONTINUE 

31  00  33  I#2, 1 1 
DX#XC-XFT*I< 

IF*DX<  32,32,33 

EVALUATION  OF  CONSTANTS  OF  INTEGRATION  FOR  CHARACTERISTICS  WHICH  ARE 
ORIGINATED  DY  THE  PROJECTILE  MOTION 

32  FRCC#-DX/<XFTSI<-XFT%I-1« 

FRC#i.-FRCC 

F2TC#F2SI-1<*FRCCGF235I<*FRC 

F6TC#F6*I-1<*FRCCEF6*I«=FRC 

TC#T*I-1<*FRCC&TSI<*FRC 

XFAC#XC/EXP*F2TC< 

GO  TO  34 

33  CONTINUE 

OBTAIN  AND  PRINT  VALUES  UF  DESIRED  OUTPUT  DATA  AT  THE  DFS I  RED 
SOLUTION  TIME 

34  F#FPRT*XFAC**VMX4 
PSI#ABS*F*SF6T-F6TC</F4T< 

THETA#PSI**VMX 1*1 2000. 

CELT  A# DOT T* THETA 
DLSTR#SFHT*THtTA 
THTT#THETA/XFTT*RNLM/12C0C.. 

X#XFAC*EXPSF2T< 

XID#X/XFTT 

XI01.-XIC 

RTHTT#XIC*THTT 

SHRT#AC**V.'1/5/<VMXl*PTMTT<**VMX3*RbFRT 

SHR#SHRT/RNIM*1000. 

VNUTL#XID*SHRT 

XCIN#XC*12. 

TC#TC* 100C . 
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X#X* 12. 

VNUSD#VNUTL*GNF 

HTC#VNUSD*VKT/DTUB/3G0. 

GD0T#HTC*^.T1T-TW< 

VNUDBfrDBF*XID**.b 

QUR#VNUSC/VNUOB 

C 

C  OUTPUT  DATA 

C  TC  TIME  ORIGIN  OF  OUTPUT  CHARACTERISTIC,  IN  MILLISEC. 

C  XC  POSITION  ORIGIN  OF-  GUTPUT  CHARACTERISTIC,  IN  INCHES 
C  HTC  LOCAL  HEAT  TRANSFER  COEFFICIENT,  IN  BTU/SEC/SQFT/DEG.R 
C  X  CURRENT  SPATIAL  LOCATION  OF  OUTPUT  CHARACTERISTIC 

C  QOOT  LOCAL  HEAT  FLUX,  IN  BTU/SEC/SQFT 

C  CUR  RATIO  OF  LOCAL  HEAT  FLUX  TO  THAT  OBTAINED  FROM  THE  DITTUS 

C  BOELTER  RELATION 

C  THETA  MOMENTUM  THICKNESS.  IN  MILS 


C  DELTA  BOUNDARY  LAYER  THICKNESS,  IN  MILS 
C  DLSTR  DISPLACEMENT  THICKNESS,  IN  MILS 
C  SHR  1000  TIMES  CF/2 

C  SHRT  EQUATION  1C2  OF  REPORT  70-18 

C  THTT  EQUATION  04  OF  REPORT  70-18 

C  RTHTT  EQUATION  104  OF  REPORT  70-18 

C  VNUTL  EQUATION  106  OF  REPORT  70-18 

C  XI  l.-X/L 

C 

WRITE*  6,8<  TC.XCIN.X, HTC, GOOT,ODR, THETA, DELTA, DLSTR 
2, SHR, SHRT, THTT, RTHTT, VNUTL, XI 
35  CONTINUE 
GO  TO  15 
60  STOP 
26  RETURN 
END 
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678. 

785. 

895. 

1020. 

1137. 

1258. 

1357. 

1450. 

1537. 

1615. 

1685. 

1757. 

1812. 

1868. 

1911. 

1951. 

1990. 

2022. 

2050. 

2077. 

2100. 

2118. 

2131. 

2148. 

2158. 

2168. 

2177. 

2182. 


2188. 

2192. 

2198. 

2200. 

2202. 

2204. 

22C6. 

2208. 

2210. 


.2  .3 

.95  .98 


.4 

.995 
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1  .0006 
1  .0007 

1  .00075 

1  .0008 
1  .0009 

1  .001 
1  .0011 
1  .00125 

1  .00145 

1  .00162 
1  .00165 

1  .0017 

1  .0018 
1  .00195 

1  .0021 

1  .00225 

l  .00243 
*STGP 
/* 


SO 


APPENDIX  B 


INPUT 

FI rst  Set  (6E10.3,  110);  1  card 

DTUB  3  Gun  tube  diameter,  inches 

ELZ  *  Chamber  length  measured  from  the  breech.  Inches 
ELD  3  Location  of  discontinuity  from  the  breech,  inches 
TW  *  Wall  temperature,  °R 

PRDL  3  Prandtl  Number  3 

CP  =  Specific  heat  at  constant  pressure,  BTU/lbm/°R 
II  a  Number  of  projectile  positions  to  be  considered 

Third  Set  (5E10.3):  II  cards 

T  3  Time,  seconds 

RO  =  Gas  density,  lbm/ft3 

T1  =  Gas  temperature,  °R 

X  =  Projectile  position  relative  to  initial  position,  inches 
V  *  Projectile  velocity,  ft/sec 

Second  Set  (3E10  =  3) ..  1  card 

VMX  *  the  exponent  m  in  Equation  52 

3  Q„2  for  1 / 7 1 h  power  law  velocity  profile 
=  0.181818  for  l/8th  power  law  velocity  profile 
3  0.166667  for  l/9th  power  law  velocity  profile 
SFG  3  Shape  factor  G  (Equation  23) 

3  1,3  for  1/7  power  velocity  profile 
AC  3  the  parameter  (eddy  viscosity)  a  in  Equation  52 
3  0.0296  for  1 / 7 th  power  law  velocity  profile 

Fourth  Set  (II,  £12.5); 

NXT  3  0  if  desired  printout  positions  have  not  been  defined 
3  1  if  desired  printout  positions  have  been  defined 
>  1  return  to  the  main  program 
TIME  3  Solution  time,  seconds 


fifth  Set*  (110).  1  card 

NN  =  Number  of  printout  positions  between  the  projectile  and 
the  breech 

Sixth  Set*  (8E10.3):  NN/8  cards 

XIC  =  Dimensionless  origin  for  characteristics  along  which 
the  output  is  desired 

OUTPUT 

First  Sequence  (II  times) 

Flrst  Set  (11  (IX  El  1,4)); 

T  *  Time*  seconds 

FI  =  Defined  in  Equation  82 

F2  =  Defined  in  Equation  81 

F4  =  Defined  in  Equation  79 

F6  *  Defined  in  Equation  122 

RNL  =  Projectile  Reynolds'  Number 

XFT  =  Projectile  position  from  breech,  ft. 

DOT  =  Defined  in  Equation  61  „  best  values  correspond  to 
curve-fitted  data  with  1/7  power  velocity  profile 

SFH  =  Defined  in  Equation  59,  but  values  correspond  to 
curve-fitted  data  with  1/7  power  velocity  profile 

SFG  =  Defined  in  Equation  60,  but  values  correspond  to  curve 
fitted  data  with  1/7  power  velocity  profile 

REFR  *  Defined  in  Equation  84 

Second  Set  (10X,  6(1 X,  El 2. 5))’ 

T  =  Time,  seconds 

RO  *  Gas  density,  lbm/ft3 

T1  *  Gas  temperature,  °R 

TFiTlfi“an3'si'xtH~ sets 'are  required  only  if  NXT  *0 
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X 


=  Projectile  position  relative  to  initial  position, 
inches 

V  -  Projectile  velocity.-  ft/sec 

HLAM  -  Projectile  base  heat  transfer  coefficient  (Equation  113) 

Se con  dSequence ; 

Fi rst  Set; 

TMS  =  Output  time.-  seconds 
XFTIN  =  Projectile  position,,  inches 

Second  Set.  Title  for  the  following  output 

Th ird  Sequence ; 

Fi fSt  Set 

TC  *  Time  origin  of  output  characteristics,  milliseconds 
XC1N  -  Position  origin  of  output  characteristics,  inches 

X  *  Current  spatial  location  of  output  characteristic 
HTC  -  Local  heat  transfer  coefficient,  BTU/sec-ft2-°R 
QDOT  =  Local  heat  flux,  BTU/sec-ft* 

QDR  =  Ratio  of  local  heat  flux  to  that  obtained  from  the 
Di ttus - Boel ter  relation 

THETA  =  Momentum  thickness,  mils 

DELTA  =  Boundary  layer  thickness,  mils 

DLSTR  =  Displacement  thickness,  mils  c 

SHR  =  Skin  friction  coefficient  (1000  times 

SHRT  =  Normalized  skm  friction  coefficient  (Equation  102) 

THTT  =  Normalized  momentum  thickness  (Equation  94) 

RTHTT  =  Normalized  momentum  thickness  Reynolds  number 
(Equation  1Q4) 

VNUTL  =  Normalized  NusseU  number  (Equation  106) 

XI  =  Normalized  longitudinal  distance  from  the  bullet 

base 
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